Drift-diffusion model of the fragmentation of the external ring structure in the 
photoluminescence pattern of indirect excitons in coupled quantum wells 
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Under optical excitation, coupled quantum wells are known to reveal fascinating features in the 
photoluminescence pattern originating from dipole orientated indirect excitons. The appearance of 
an external ring has been attributed to macroscopic charge separation in the quantum well plane. 
We present a classical model of non-linear diffusion to account for the observed fragmentation of 
the external ring into a periodic array of islands. The model incorporates the Coulomb interactions 
between electrons, holes and indirect excitons. At low temperatures, these interactions lead to pattern 
formation similar to the experimentally observed ring fragmentation. The fragmentation is found to 
persist to temperatures above the quantum degeneracy temperature of indirect excitons. 
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Since the first theoretical predictions that excitons - 
bound states of electrons and holes - may undergo Bose 
Einstein condensation (BEC) [1], intensive studies have 
sought to realize the phenomenon experimentally. This 
arduous journey has inspired a diverse array of exper- 
iments and investigations GHSl. Novel structures and 
techniques have been employed which, regardless of 
their success or failure, have deepened our understand- 
ing of the excitonic properties of semiconductors. 

A system which remains a strong candidate for ob- 
serving exciton BEC is that of indirect excitons in cou- 
pled quantum wells (CQWs) [2,3]. The samples consist 
of two thin layers of semiconductor material stacked 
between layers with a higher band gap energy. A volt- 
age applied in the growth direction shifts the electron 
and hole energies with respect to each other. In this 
arrangement, optical excitations create indirect excitons 
where electrons and holes are confined to adjacent quan- 
tum wells by tunneling through the barrier region. The 
spatial separation of electron and hole wave functions 
extends the optical lifetimes of excitons well beyond the 
time needed for thermalization to the lattice temperature 
HTHTTl. These structures provide optimal conditions for 
creating a cold dense exciton gas - a prerequisite for BEC. 

Some of the first results from this type of experiment, 
published nearly a decade ago 0131, revealed some re- 
markable features in the spatial photoluminescence (PL) 
pattern of indirect excitons. These include the inner ring 
and localized bright spots. Of interest here is the bright 
external ring which can be greater than 100 /im in diame- 
ter and is separated from the central bright spot around 
the laser by a predominantly dark region. Previous stud- 
ies IH2HT4I have shown that the ring forms at the inter- 
face between spatially separated electron-rich and hole- 
rich regions. This is due to a difference in the capture 
efficiency of the injected electrons and holes which re- 
sults from their different masses. An abundance of holes 
builds around the laser spot as some fraction of electrons 
leaks to the electrodes. In the presence of an electric field, 



the doping in the electrode layers leads to a background 
electron concentration within the wells. The result is 
a pool of holes surrounded by a sea of electrons with 
excitons forming at the boundary. 

Perhaps the most remarkable feature observed is the 
fragmentation of the external ring into a periodic array of 
islands for temperatures below « 3 K |3j[l5Hl6|- We note 
that a previous model [17] explains the macroscopic or- 
dering via a process of stimulated scattering of excitons 
into the ground state whilst another model [18] assumes 
a BEC present in the ring. Both imply that fragmentation 
is a signature of quantum degeneracy. 

An alternative model explained the ring fragmenta- 
tion in terms of attractive interactions between excitons 
lH9l l20l . In this framework, a ring of uniform density 
becomes unstable and the attractive interaction leads to 
droplet formation. However, measurements of the blue 
shift in the exciton emission [15] reveal clear evidence 
that the interaction is repulsive and that this mechanism 
cannot be responsible for the effect. 

In this letter, we present a model to asses the impor- 
tance of Coulomb interactions in the external ring. We 
find that at low temperatures, these interactions lead 
to periodic modulation of the exciton density. Whilst 
spatial coherence measurements suggest a statistically 
degenerate exciton gas is present in the ring HUE], we 
illustrate that fragmentation can occur due to classical 
mechanisms and cannot unambiguously be taken as ev- 
idence of degeneracy. Studies of the role of Coulomb 
interactions in this system have been attempted previ- 
ously l22l 1231 and have highlighted their effect on the 
appearance of the external ring. However, due to the 
computational complexity of the problem, these works 
have been limited to the case of ID geometry or to a low 
density regime and are insufficient to capture fragmen- 
tation of the ring. In Ref. |23 |, it is shown that inclusion 
of Coulomb terms is necessary to explain the threshold 
dependence of the ring radius on excitation power and 
confirms the significance of the contribution made by the 
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interactions included here. 

Conceptually, we find that at low temperatures, the 
diffusive nature of carrier transport diminishes. In this 
regime, carriers move under the action of strong electric 
fields provided mainly by the in-plane charge separa- 
tion of electrons and holes. Under these conditions, a 
ring of uniform density is an unstable configuration and 
the system naturally evolves to a fragmented ring as a 
process of energy minimization. 

Based on the approach used in mi EJ [23j EU, the 
following set of coupled equations were used to model 
the creation, transport and decay of carriers in CQWs: 



~dt 

dt 

dn x 

dt 



= V [D e Vn e + n e n e VU e ] + 



nl 0) - ft e 



+ A e - K7ft e fth, 



— = V [D h Vft h + ^hfthVLTh] + A - wn e n h , 



= V [D x Vft x + ^ x ft x VU x ] + wn e ti h - 



(1) 



The solution is the density distributions ft z of electrons, 
holes and excitons (i = e,h,x respectively). The first and 
second terms inside each of the square brackets are the 
diffusion and drift currents respectively. D z and \i{ are 
the diffusion coefficients and mobilities. Each diffusion 
coefficient is assumed constant and we use the classical 
limit for the mobilities, \i\ « D z 7(fcgX). Far from the laser 
induced heating at the excitation spot, the temperature 
of each species is defined by the lattice temperature, T. 
Coulomb interactions are included via the potentials U e , 
Uh and U x . For electrons and holes residing in adjacent 
quantum wells, one has 

li e = ft e *V - n h *V d + n x *(V -V d ), 

U h = n h *V - n e *V d + n x *(V -V d ), (2) 

U x = (ft e + ft h + 2ft x ) * (V - V d ), 



where V a (r) = e 2 /(4:Ti£o£ r ^/|r| 2 + a 2 ) and V denotes the 
convolution; f * g = J f(x - t')g{t')dr f . Here, r is 
the in-plane coordinate. All calculations are made for 
GaAs/Al .33Ga .67As CQWs studied in Ref.|3|. We use 
d = 12 nm for the separation between electron and hole 
layers and dielectric constant e r = 12.9. Both d and e r are 
critical parameters in the onset of fragmentation since 
they scale the strength of the interaction terms. Macro- 
scopic charge separation creates an in-plane potential 
gradient under which electrons and holes drift towards 
the ring position. The separation of opposite charges into 
adjacent wells leads to repulsion in the exciton-exciton 
interaction and also in the exciton-electron and exciton- 
hole interactions. The interactions make a significant 
contribution at low temperatures where transport due to 
Coulomb forces dominates over diffusive mechanisms. 

The remaining terms on the right hand side of Eq. (1) 
account for the creation and decay of each species. The 
homogeneous source term for electrons, which is due to 
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FIG. 1: Spatial density distributions for electrons (a), holes 
(b) and indirect excitons (c). An optical pump of FWHM 
10 f/m is focused at the origin with peak generation rate 
of 2 x 10 10 cm _2 ns _1 . Background electron density nf* = 
0.5 x 10 9 cm -2 . Other parameters are given in the main text. 
A cut along the x-axis of the density distributions (d) and po- 
tential (e). The pump profile is shown by black dots. 



a constant flux through the CQWs, acts to restore a back- 
ground density, nf } in a characteristic time T e . Both n\ ' 
and x e depend on the applied electric field. To reduce 
the number of control parameters, we modeled control 
of the electric field by varying nf^ and fixed T e = 50 ns 
which is close to the experimentally determined value 
in Il25l . A(r) is the injection rate of unbound electron- 
hole pairs into the CQWs by the laser. The imbalance in 
the generation rate of electrons and holes can be incor- 
porated by A e < A. Without loss of generality, we set 
A e = and compensate by reducing A. This removes the 
central bright spot from simulations. The binding rate 
of free pairs into excitons is proportional to the overlap 
of the electron and hole densities. The fit parameter w 
is inversely proportional to the exciton formation time 
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and fragmentation is sensitive to its value. Here we used 
w = 10 3 cm 2 s~ 1 l26l . This leads to an exciton density of 
about 10 10 cm~ 2 which has been determined experimen- 
tally via the exciton blue shift [3|. Reducing w, reduces 
n x and diminishes the interactions, preventing ring frag- 
mentation. This term is a decay channel for electrons and 
holes and a source term for excitons. t x is the optical life- 
time of excitons and is nearly constant with respect to n x 
IITTl . To simplify the model, we fixed t x = 50 ns and have 
not explored the effect of its dependence on temperature 
or electric field. D x = 0.2 cm 2 s _1 was used for the exci- 
ton diffusion constant. This is consistent with the model 
used in |TT) where the effect of a QW disorder potential 
of « 1 meV was included. The electron and hole diffu- 
sion constants are D e = 30cm 2 s _1 and Dh = 15cm 2 s _1 , 
comparable to experimentally measured values in [25 J. 

Figs, la-d show the density distributions satisfying (1- 
2). The equations were solved dynamically with the po- 
tentials Mi being continuously recalculated on each time 
step until a steady state solution reached. The simula- 
tions use initial conditions of uniform density n e = nf^ 
and fth = = 0. At a circular boundary well beyond the 
ring position, the electron density is fixed to nf* and the 
hole and exciton densities to zero. A Gaussian excitation 
profile for A(r) is focused on the center of the sample. The 
results show a periodic modulation of the exciton den- 
sity along the ring. Additionally, a slight modulation of 
the ring radius is seen which is not observed in exper- 
iments. Figs. Id and le show cross sections along the 
line y = for the density distributions and potentials re- 
spectively. In-plane charge separation causes an E-field 
of « 2eVcm _1 driving electrons and holes to the ring 
position thus enhancing the generation rate of excitons. 

A feature of the experimental data which is well repro- 
duced by the model is the dependence of the ring frag- 
mentation on the lattice temperature. Fragmentation is 
observed only below a critical temperature |3, 15|. This 
has been interpreted as the indirect exciton degeneracy 
temperature where the statistics cross from the classical 
to the quantum regime [17 \. In our model, the lattice tem- 
perature T appears in the denominator of the interaction 
terms which, as a result, become greater in magnitude 
than the diffusive terms for V17 Z > fcBT(Vn z )M- 

The inset in Fig 2a shows the exciton density profile at 
the ring position, r ring for various temperatures. The ring 
position is defined as the location of maximum exciton 
density for each angle 6 about the excitation spot. Ring 
fragmentation appears with decreasing temperature as 
the drift currents dominate over the diffusive currents. 
This is illustrated in the main panel where the average 
contrast between the peaks and dips along the ring is 
plotted against T. The model captures the temperature 
dependence of the pattern formation observed in exper- 
iments l3lH5l. The onset of fragmentation is abrupt and 
a critical temperature can be identified. This temper a- 
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FIG. 2: a) Ring contrast, <(J™ ax - I™ n )/Ip L ax > against temperature 
for various nf\ Average ring radius is kept constant by ad- 
justing A. Inset: PL intensity profiles J PL along the ring for the 
points on the red curve shown by arrows, b) Average density in 
the external ring against temperature. The dashed line marks 
the transition from a uniform to a fragmented ring. The solid 
line marks the crossing from classical to quantum statistics. 

ture is deceivingly low and could be misinterpreted as 
the degeneracy temperature. The critical temperature 
varies with exciton density which can be controlled by 
adjusting A and simultaneously whilst maintaining 
a fixed average radius. Each curve corresponds to a dif- 
ferent value of nf* with A chosen to give a radius of 
50 fim. The seemingly noisy data in Fig 2a is a conse- 
quence of the state of the ring evolving to configurations 
with different numbers of islands. The degenerate so- 
lutions have slightly different contrasts and the effect is 
most pronounced at the highest density. 

Fig 2b shows the spatially averaged density in the ring 
against temperature for different values of nf\ The 
dashed line marks the average density at each critical 
temperature. Also shown is the line T = Tq. Here, 
T = (nh 2 n x )/(2M x kB) is the exciton degeneracy temper- 
ature (M x = 0.22rao is the exciton mass). The model 
predicts that fragmentation can occur both well above 
and below Tq. At the highest densities examined, an in- 
crease in exciton density is observed as the temperature 
is further reduced beyond the onset of fragmentation. 
In PL experiments, this manifests itself as an increase in 
the blue shift in the exciton spectrum which has been ob- 
served in Ref . 1 15 ]. In that work, the blue shift was found 
to also increase with temperature when above the critical 
temperature. Our model would be consistent with this if 
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FIG. 3: Exciton density distribution due to a line excitation 
for A rand = (a) and A rand = 1 jueV (b) at T = 1 K and nf } = 
10 9 cm -2 . The density along the external line with and without 
disorder at T = 1 K (c). Periodic boundaries are at y = ±128/im 
and n£ 0) = 10 9 cm- 2 . 



the temperature dependent t x described in |9] was used. 

Due to doping, there is an abundance of free electrons 
in the layers surrounding the quantum wells. Arguably, 
these electrons may accumulate around the excitation 
spot to screen the in-plane potentials U e and which 
would suggest that Coulomb interactions are not promi- 
nent enough to be responsible for the fragmentation of 
the ring. However, on comparison of available experi- 
mental data, an alternative picture presents itself. Ring 
fragmentation is seen in the work of Butov et al. [3], but 
not in the work by Snoke et al\2\, where it is estimated 
that a much greater density of free electrons is present 
in the layers adjacent to the CQWs. From this, we can 
deduce that in the former there are insufficient free elec- 
trons to screen the Coulomb interactions and prevent 
ring fragmentation. The detailed effect of a free electron 
gas in the adjacent layers is beyond the scope of this work 
and its effect on pattern formation is an open question. 

To further understand the physical mechanism that 
leads to ring fragmentation, it is instructive to exam- 
ine the case of a line excitation. In this geometry the 
laser is focused not to a single point but to a narrow 
line extended across the entire sample. This leads to the 
appearance of two parallel external lines in the indirect 
exciton PL pattern, either side of the laser. In Fig. 3 we 
present results for this situation. The laser has a Gaus- 
sian profile in the y direction and is homogeneous in the 
x direction. Periodic boundary conditions are used to 
simulate an infinitely extended line. Unlike the external 
ring in the point excitation geometry, the external line 
does not spontaneously fragment. Instead, a small per- 
turbation is needed to observe the effect. This is included 
by adding a random disorder potential to either I7 e , 
or U x . The density distributions in Figs. 3b and 3a are 
from simulations with and without a perturbation via 
the electron potential U e , respectively. The amplitude of 
the disorder potential used, A ran d = 1 f/eV is orders of 
magnitude smaller than the fluctuations in U e along the 



external line and the correlation length of the disorder 
(2/im) is a few times less than the period of density mod- 
ulation. This confirms that the density modulation is 
not correlated to the disorder but disorder is required to 
evolve the state of the line from the metastable solution 
shown in Fig 3a to the stable solution shown in Fig 3b. 
In the circular geometry, the numerical discretization of 
Eqs. (1-2) onto a rectangular grid provides some implicit 
distortion that triggers the fragmentation. 

Fig. 3c shows the exciton density at the line position 
with and without disorder. A consequence of fragmen- 
tation is the reduction in the exciton density and, there- 
fore, a lowering of the energy density associated with 
the dipole repulsion between excitons. We can conclude 
that in the simulations, the external line buckles into a 
wavy line in order to redistribute excitons over a larger 
area and reduce the energy density. We found that any 
non-zero perturbation triggers an energy minimization 
process, leading to fragmentation of the line. 

The model presented provides insight into the role 
of Coulomb interactions in the formation of the exter- 
nal ring in the indirect exciton PL pattern. We have 
demonstrated that within a classical framework, peri- 
odic modulation of the exciton density on macroscopic 
length scales can occur. In the classical picture, fragmen- 
tation is a process of redistributing charges to minimize 
the potential energy associated with dipole repulsion. 
We found strong similarities with the available experi- 
mental data. In particular, the calculated temperature 
dependence of the ring contrast is in strong agreement. 
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